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ABSTRACT 

Chandra's high resolution observations of radio galaxies require a revisit of 
the relevant electron acceleration processes. Although the diffusive shock par- 
ticle acceleration model may explain spectra of spatially unresolved sources, it 
encounters difficulties in explaining the structure and spectral properties of re- 
cently discovered Chandra X-ray features in several low-power radio sources. We 
argue that these observations strongly suggest stochastic electron acceleration 
by magnetized turbulence, and show that the simplest stochastic particle accel- 
eration model with energy independent acceleration and escape timescales can 
overcome most of these difficulties. We use the bright core of the western hotspot 
of Pictor A as an example to demonstrate the model characteristics, which may 
be tested with high energy observations. 

Subject headings: acceleration of particles — galaxies: individual (Pictor A) - 
galaxies: jets — radiation mechanisms: non-thermal 



1. Introduction 

In the classical diffusive shock particle acceleration model the observed relativistic elec- 
trons in radio galaxies are accelerated at the edges and knots of jets. For relatively more 
luminous FR II radio galaxies, electrons may also be energized in radio hotspots, where jets 
terminate. It is usually assumed that the jet power is initially stored in the bulk motion of 
outflows and is conve rted into the internal ene rgy of the jet plasma by shocks at the observed 



radio bright features (jBlandford fc Reedll974l ). The model requires that particles repeatedly 



pass the shock front, characterized by a jump in the density and velocity field, to reach very 
high energies. Although the physical processes responsible for the scattering and/or diffusion 
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of particles through the shock front in most collisionless astrophysical plasmas are not well 
established, it has been widely cited due to it s simplicity and its achi evement in explaining 



the spectra of spatial unresolved sources (e.g. lHardcastle et al. Il200ll ). 



It however is challenged by recent high spatial and spectral resolution observations with 
Chandra. Although nonthermal X-rays from radio galaxies can be produced by relativistic 
electrons through synchrotron, synchrotron self-Comptonization (SSC), and inverse Comp- 
ton scattering of the cosmic microwave background radiation (CMBR) or other background 
photons, X-rays from jets of FR I galaxies and s ome low-power hotspots of FR II galaxies 



are generally attributed to synchrotron process (jWorrall. Birkinshaw. fc Hardcastle 1 12001 



Hardcastle et al. 1 12004 ) . If the magnetic field is nearly in energy equipartition with the 



relativistic electrons, the synchrotron cooling timescale of electrons producing the observed 
X-rays is short, and the electrons need to be accelerated very close to the observed compact 
X-ray features, which usually are displaced from peaks of the corr esponding radio sources 



and occasionally do not have optical and/or radio counterparts (lHardcastle et al. 1 12001 



2007al ; Marshall et al. 1 120021 ; IWilson fc Young Il2002l ). The diffuse X-ray emission from jets 
of FR I galaxies, on the other hand, indicates a distributed particle acceleration process. The 
X-ray spectral hardening and relatively high X-ray flux as compared to the extrapolation 
from optical fluxes of some compact sources are also puzzling. 

The bright core of the western hotspot of Pictor A (WHPA), a famous nearbdj] FR II 



galaxy , is an excellent candidate to demonstrate diffic ulties of the shock model. iMeisenheimer. Yates fc Roe: 
( 119971 ) discussed its radio to optical spectrum, and IWilson. Young fc Shopbell I (120011 ) re- 
ported X-ray observations with Chandra. The latter studied in detail several possible X-ray 
emission processes and found that none of them can explain the broadband spectrum (BBS) 
naturally. They concluded that the X-rays are most likely synchrotron emission by an elec- 
tron population distinct from these responsible for the radio to optical spectrum. In this 
picture the fact that the radio, optical, and X-ray emissions come from the same spatial 
region would have to be a coincidence. 

In this Letter, we show how the simplest stochastic particle acceleration (SA) model 
explains the BBS of WHPA. The model also predicts a spectral cutoff near the Chandra 
band and an intrinsic volume of 0'.'03 3 of the particle acceleration region (AR), which may be 
tested with future observations. In § 2, we discuss the major difficulties of the shock model 
in explaining Chandra observations and why the SA by turbulence is favorable. A simple 



1 At a redshift of z = 0.035, the luminosity distance to Pictor A is c?l = 151.9 Mpc, and the conversion 
scale is 1" = 0.687 kp c for the modern cosmology with Hm = 0.27, £!a = 0.73, and Hq = 71kms _1 Mpc -1 
jSpergel et al. II2003I ). 



- 3- 



SA model is built and applied to WHPA in § 3. In § 4 we draw conclusions and discuss how 
the model may also explain other Chandra X-ray features. 



Shock VS Stochastic Acceleration 



The synchrotron cooling time of an isotropic electron population is given by r syn = 
9mgC 5 /4e 4 .B 2 7 = 24.G'-f~ 1 B~ 2 yr, where m e , e, 7, c, and B are the electron mass, charge, 
Lorentz factor, the speed of light, and magnetic field, respectively. Here and in what follows, 
c.g.s. units are adopted. The synchrotron power spectrum peaks near ~ 1.7 x 10~ n B~f 2 
keV. For a 0.1 mG magnetic field, typical for jets and hotspots of radio galaxies, the X- 
ray emitting electrons have 7 > 3 x 10 7 and r syn < 82 yr. The corresponding light travel 
distance is < 25 pc, which is comparable to or smaller than the Chandra resolution for 
most radio galaxies. The X-ray emitting electrons therefore need to be accelerated near 
the observed X-ray sources. In the standard diffusive shock model, particles are accelerated 
at a thin shock front and are carried away by the shocked flow into the downstream. The 
simplest 1-D model predicts a power-law synchrotron spectrum that cuts off at the frequency 
where the synchrotron cooling time is equal to the flow travel time from the shock front 
( iHeavens fc Meisenheimed 119871 ). Then the spatially integrated spectrum from any segment 
of the downstream region can be fitted with a broken power-law, which cuts off at the cutoff 
frequency of the boundary closer to the shock front. The break frequency corresponds to 
the cutoff frequency of the farther boundary. The model has difficulties in accounting for 
the following Chandra observations: 

(i) If the diffuse X-ray emission from jets is to be associated with spatially unresolved 
shocks, the jet dynamics must be sophisticated enough to induce many shocks, (ii) The 
X-ray spectral hardening is observed in a few jet knots of FR I galaxies and hotspots of FR 
II sources, e.g., knot WK7.8 of PKS 0637-752 JSchwartz et al. Iboooh. some knots in the 



lets of M87, 3C 120 and NG C 6251 flWilson fc Young 



2002 



Harris. Mossman. &: Walker 



20041 ; ISambruna et al. 1120041) WHPA (I Wilson. Young fc Shopbell ||200jJ) , and hotspot Pl/2 
of 3C227 JHardcastle et al. I l2007ah . bermer fc Atovan I J2OO2F suggested that it may be 
due to the Klein-Nishina effects of inverse Compton scattering of CMBR by relativistic 
electrons that produce the observed X-ray through synchrotron process. The model requires 
an extremely low magnetic field with an energy density comparable to that of the local 
CMBR or very strong Doppler effects and may not be applicable to most of the observed X- 
ray sources in radio galaxies, which are still quite powerful and do not show strong Doppler 
effects. The extreme hardening of hotspot Pl/2 of 3C227 and inner (< 50 / . / ) jet of Centaurus 
A is even more challenging. The photon spectral indexes of these sources can be as small as 
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1.6 corresponding to an electron spectral index of 2.2, which can not be a cooling spectrum in 
the shock model. The lifetime of these sources therefore must be shorter than the synchrotron 
cooling time of X-ray emitting electrons. Given the observed source extension, a thin shock 
front may not be able to accelerate enough electrons to produce the observed X-ray flux 
over such a short time scale (< 100 yr). (iii) The displacement of X-ray peaks with respect 



to th e corresponding optical and/or radio peaks (IWilson fc Young 1 120021 ; lHardcastle et al. 



200ll ). The shock model predicts that higher energy emission is produced closer to the 
shock front, and the distance of synchrotron emission peak from the shock front is inversely 
proportional to the square root of the emission frequency for a constant downstream velocity 
and magnetic field (IMeisenheimer et al. Ill989l ). There is no observational evidence for such 
a source structure. On the contrary, much more complicated structures are seen in a few 
low-power hotspots of FR II galaxies (lHardcastle et al. Il2007al ). 



While radio observati ons are generally consistent with shock models for hotspots (IMeisenheimer et al 



19891 ; ICarilli et al. I Il99ll ). there are in dications of electrons being accelerated directly by 
magnetic fields in giant radio galaxies (IKronberg et al. II2004I ). Optical detections of some 
of these radio sources sho w that electrons are accelerated in ext ended regions that may not 
be associated with shocks (IMeisenheimer. Yates fc Roeser 1119971 ). Little is known about this 
distributed acceleration process driven by free energy dissipations. The high Reynolds num- 
bers of the observed radio sources suggest that turbulence may mediate these dissipations 
and accelerate some particles, the so-called SA, a second order Fermi acceleration that has 
been ignored in most astrophysical situations due to its relatively lower efficiency The ob- 
served cutoffs of synchrotron spectra suggest that t he acceleration i s less efficient than that 



given by the diffusive shock with Bohm diffusion (IBrunetti et al. 1 120031 ; lHardcastle et al. 



20041 ). Moreover, for strongly magnetized plasmas, t he Alfven velocity can be comparable 
to c and SA can be effective (jPetrosian fc Liu Il2004j ). These combined with the Chandra 
observations strongly suggest that SA may be the dominant acceleration mechanism in radio 
galaxies. 

SA occurs wherever there are turbulent magnetic fields and can operate over an extended 
region with its size dictated by the turbulence generation and decay rates. Emission from the 
AR therefore can be significant. We consider the simplest SA model with energy independent 



acceleration (r ac ) and escape (r ( 



times. Then the steady-state particle distribution in th e 



AR is given by (jSchlickeiser I Il984j ; iPark fc Petrosian I Il995l ; iLiu. Melia fc Petrosian 1 120061 ) : 



iV ac (7) oc 7 d+2 exp(-7/ 7c ) 

J °° x 5 ' 1 (1 + xf +s e^ x '^dx for 7 > 7 inj 
Jq 1 x 6 - 1 ( 1 - x) 3+5 e 7X / 7c dx for 7 < 7 inj 
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where 5 = (9/4 + Tac/resc) 1 / 2 — 1.5, y c = 9rrtgC 5 /4e 4 i? 2 r ac , and we have assumed that electrons 
are injected into the AR at 7i n j and the synchrotron cooling dominates. 

As in the shock model, electrons escaping the AR may produce most of the observed 
emission in a region dominated by cooling processe s. The corresponding particl e distribution 



at time t since the acceleration starts is given by (IBfumenthal fc Gould Ill970l ) 

AT f* 7(7o) Wac(7o) 

Jo 717) Tesc 

where 7(7) < is the energy loss rate and t — t = drf/ytf). Note that if t — t < 
dy'/y^'), the integrant should be set to zero. The thin dashed and dotted-dashed lines 
in Figure [TJ give iV ac and iV cool , respectively, for 7 in j = 1800, 5 = 1.48, y c = 1.04 x 10 7 , 
t = 17.9r ac (see § |3]on values of these parameters) and 7 oc — -y 2 . We see that the spatially 
integrated particle distribution as indicated by the thin solid line is very similar to the result 
of the shock model, although the AR can be much more extended in the SA model. 



3. Application to WHPA 



To apply the model to ho tspots of radio galaxies, one needs to tak e into account the 



effects of adiabatic expansion (jCarilli et al. I Il99ll ; iBrunetti et al. II2003I ). Such an expan- 



sion is expected since the pressures of hotspot plasmas are usually much higher than their 
surroundings. To simplify the model, we assume that this expansion is decoupled from the 
acceleration and cooling processes, i.e. the AR experiences a quick adiabatic expansion be- 
fore merging into the cooling region (CR). For a spherical AR with its radius increased by 
a factor of Ar during the expansion, the magnetic field energy density decreases by a factor 
of Ar 4 and the energy of relativistic electrons decreases by a factor of Ar. After the expan- 
sion, the plasma enters a relatively uniform cooling region. We are mostly interested in the 
synchrotron emission. To reproduce the optically thin synchrotron spectrum of the AR with 
the magnetic field of the CR, one needs to increase the electron energy and total number 
by a factor of Ar and Ar 2 , respectively, since the synchrotron frequency and luminosity are 
proportional to -B7 2 and -B 2 7 2 , respectively. The thick lines in Figure [TJ show these equiva- 
lent electron distributions for Ar = 2.5. Emission from the adiabatic expansion phase can 
be ignored as far as our assumption of quick expansion withheld. We see that the overall 
spectrum becomes harder right below the cutoff due to dominance of the AR there. 



Figure [2] shows the BBS of WHPA, where all data are taken from lWilson. Young fc Shopbell 

(120011 ) . To fit it with the synchrotron spectrum of the electron distribution in Figure [TJ 



one can have a Ar > 1 and use the hardening of the electron distribution below the cut- 
off to produce the observed X-ray spectral hardening. For a given B, one can constrain 
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7c and Ar by fitting the X-ray spectrum. The radio spectral index a = 0.74 implies 
5 = 1.48. Then the frequency of the synchrotron spectral break from the radio to X-ray 
ban d determines t. The intrinsic dime nsion of the bright core of WHPA is 2'! x 1". x 



1" (jPerley. Roser &: Meisenheimer 1 119971 ). the equivalent radius R = 0'.'78 = 1.65 x 10 21 



cm. The normalization of the elec tron distribution is th en determined by the flux density 



(Doppler effects are ignored here). iMigliori et al. I (120071 ) find the electron to magnetic field 



energy density ratio U e /UB ~ 25(7 in j/50)~ a6 in the radio lobes of Pictor A. Since relativistic 
electrons and magnetic field experience the same adiabatic expansion loss from hotspots to 
radio lobes and the former is also subject to radiative losses, this ratio should be higher in the 
hotspot. We assume U e /UB = 15(7i n j/1800) 1_<5 and adjust B to fit the BBS. The thick solid 
line in Figure [2] shows the corresponding best fit, where, besides the parameters discussed 
above, B =85 fiG and the electron density n e = 2.8 x 10~ 6 cm -3 in the CR. The X-ray 
spectral hardening is reproduced with a Ar = 2.5. The lifetime of the CR t = 5.2 kyr. The 
energy flux carried away by electrons and magnetic field from the AR is 1.5 x 10 45 erg s^ 1 , 
which is a reasonable low limit of the jet power. For an electron-proton plasma without a low 
energy component, the corresponding Alfven velocity is 0.26c implying a highly magnetized 
plasma. 

The equivalent radius of the AR R a = 4.0 x 10 19 cm. The corresponding light crossing 
time of 42 yrs is much longer than the cooling time of 7.4 yr at 7 C (i.e. the acceleration 
time r ac ) in the AR. Therefore the shock model has difficulties in reproducing the X-ray flux 
since the lifetime of the X-ray emitting electrons (with 7 ~ j c ) is too short for a thin shock 
front to generate enough such electrons to reproduce the observed flux. One can remedy the 
shock model by decreasing the magnetic field in the downstream region or increasing the 
shock front cross section. The former requires a magnetic field energy density far below the 
equipartitional value with the electrons, which seems unlikely given that a strong magnetic 
field is required to convert the jet energy into the energy of particles. For a downstream flow 
velocity va, the latter requires a cross section greater than 47ri2^/3u c jT ac ~ (0.1c/fd)(0'/3) 2 . 
From 5 = 1.48, we have r esc = r ac /6.6, the thickness of the AR needs to be on the order of 
1.1 lyr ~ 0.5 marcsecond. The corresponding cross section is 2.5 x 10 41 cm 2 ~ 0'.'24 2 , which 
is smaller than that of the shock model. 

The adiabatic expansion phase plays a crucial role in the X-ray spectral hardening. The 
synchrotron spectra for Ar=1.0, 2.0, and 3.0, are shown by the thin dotted, dotted-dashed, 
and dashed lines in Figure [21 respectively, where n e is proportional to Ar s due to the shift of 
the injection energy by the adiabatic expansion. The other model parameters are the same 
as the fiducial model. We see the spectral hardening is prominent only for Ar > 2 and there 
is no hardening without adiabatic expansion, i.e. Ar = 1. Because 7; n j is fixed, the radio 
spectrum also changes slightly with Ar. Electron injection from jet to a hotspot AR has long 
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7inj = 100 bv ICarilli et al. I (jl99l|), and 7 inj = 1000 bv lWardle et al. I (Il998f l. In the context 



been debated, the so-called injection problem of shock models ( e.g.. ICarilli et al. Ill99lf). A 
wide range o f is used by diffe rent authors, e.g. 7i nj = 10 used by ICroston et all (120051 ) 



of SA, an electron energy comparable to the proton rest mass ene rgy is favored bec ause lower 
energy electrons can be accelerated efficiently by whistler waves (ILiu et al. 1 120061 ) . A recent 
reanaly sis of the BBS of the h otspots in Cygnus A also suggests a sharp spectral break at this 
energy (jStawarz et al. 1120071 ). We therefore choose 7 in j = 1800 for the fiducial model. The 
two thin solid lines in Figure [2] show the spectra for 7 in j = 10, and 3000. The latter under- 
predicts the 90 cm radio flux by ~ 60%. Here n e oc and U e /UB oc 7 i 1 n 7 <5 '. Figure [2] also 
shows the spectral evolution assuming a constant injection power and constant magnetic field 
and electron density in the CR. As expected, the total luminosity of the hotspot increases 
with time. At t — 0.4 Myr, the X-ray flux is dominated by the SSC component, which may 
expla in the observation that the X-ray emission of high-power hotspots are dominated by 
SSC (IHardcastle et al. 1 120041 ) . The SSC component becomes more prominent if a bigger 
Ue/Us is chosen for the fitting, which also results in a lower B, higher j c , and higher energy 
flux from the AR. 

Several models have been proposed for the X-ray spectral hardening before. None of 
them appears to be applicable to WHPA. We also consider the Klein-Nishina effects of SSC 
on the electron distribution. Figure 2 shows that the volume of the emission region needs to 
be at least 10 times smaller than the fiducial model to make the SSC coolingdominant. This 
leads to an electron to magnetic field energy density ratio greater than 150 □ and significant 
SSC contributi ons to X-ray emiss ion. The SA itself can also produce hardening near the high 
energy cutoff (jSchlickeiser II19841 ). However, this is significant only for r csc ^> r ac (5 ~ 0). 
The corresponding electron spectral index is ~ 1, which is not consistent with observations. 



4. Conclusions and Discussion 

We have shown that SA reproduces the X-ray spectral hardening of the bright core of 
WHPA and predicts a spectral cutoff near 1 keV. With few modification, the model can 



2 The synchrotron luminosity L syn oc NqR 3 B 2 ^ x , where Nq is the normalization of the electron distribu- 
tion and jx is the Lorentz factor of electrons producing keV X-rays through synchrotron process, i.e. B^\ 
is on the order of 10 11 G. To make the SSC cooling dominant, the energy density of synchrotron photons 
needs to be greater than the magnetic field energy density, i.e. B 2 R 2 < 2L syn /c. The electron energy density 
is proportional to iV o7i 1 n T 5 . Then we have B 2 /N ol ^ 5 oc B 4 R 3 j%/L syn ^ s < 10 n (2/ 'cfl 2 Lj^ygr 1 G. One 
needs to increase 7i n j, which is constrained by radio flux densities, to bring the magnetic field and electrons 
close to energy equipartition. 
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be applied to other Chandra X-ray features with similar properties. We emphasize that 
the shock model has difficulties in explaining these observations since the lifetime of X-ray 
emitting electrons is too short for the shock to produce the observed high flux levels. This 
challenges any models with particles accelerated efficiently in very small regions With a very 
large shock cross section (> 0'/3 2 ) and high downstream flow velocity (> 0.1c), the shock 
model may fit the BBS of WHPA by assuming a quick adiabatic expansion phase of the 
downstream flow, similar to the SA model discussed above. However, a power-law spectrum 
with a photon spectral index of 2.24 is expected beyond the X-ray band except that the 
acceleration rate of highest energy electrons is fine-tuned to be equal to the cooling rate of 
X-ray emitting electrons. In this case the model predicted BBS will be indistinguishable from 
that of the SA model. For Bohm diffusion , the acceleration rate is equ al to eBAr 2 c/27r, the 



emission spectrum cuts off at ~ 70 MeV (IBlumenthal &: Gould Ill970l ). Observations above 



the Chandra band (e.g. with Suzaku) can readily distinguish these two models. 

The SA is less efficient and operates over an extended region with its size determined 
by the turbulence generation and decay processes. We predict that the observed harder X- 
ray emission originates from the AR while the radio (optical) emission comes from the CR. 
The energy dissipation and the accompanying SA may still be triggered by shock waves. In 
this case, the traditional study of shock dynamics can be used to make predictions on the 
displacements between peaks in the X-ray and radio (optical) images. These may explain 
observations of bright knots in FR I galaxies. The dynamics of jet flows is sophisticated. 
Some of the energy dissipation regions may not be associated with shock fronts. In this case, 
one may use MHD simulations to identify regions with strong energy dissipation and apply 
the SA model to make predictions on the source morphology. MHD simulations may also 
address the dynamics related to the turbulence and the adiabatic expansion phase introduced 
to produce the observed spectral hardening. 
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Fig. 1. — Effective electron distributions for Ar = 2.5 (thick) and 1.0 (thin). The dashed 
and dotted-dashed lines are for the acceleration and CR, respectively. The solid lines are 
the sums of the two. The normalization is for the fiducial model. See text for the rest of the 
model parameters. 
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Fig. 2. — Best fit to the spectrum of W HPA and effects of Ar, 7^, and t. The data are 
from IWilson. Young fc Shopbell I (120011 ). The thick solid line is the best fit with Ar = 2.5, 
B = 85 /iG, n e = 2.8 x 10~ 6 cm~ 3 , and t=5.2 kyr. The thin dashed, dotted-dashed, dotted 
lines are for Ar = 3.0,2.0,1.0, respectively. The two thin solid lines are for 7i n j=3000, 
and 10 with the former having a lo wer radio flux. The downward t riangl e radio points are 
for spatially unresolved fluxes from iPerley. Roser &: Meisenheimer I (119971 ). Their relatively 
harder spectrum suggests a spectral break near 7 ~ 1000. The higher energy components 
are produced through SSC. The thick dotted-dashed and dashed lines are for i=2.3 and 400 
kyr, respectively. The SSC component dominates the X-ray flux for the latter. 



